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Abstract 


A straight, single channel membrane humidifier was constructed to measure temperature and moisture profiles along both the donor and receiver 
channels. A persulfonic Nafion membrane was used as the water exchange medium. 

We report on results obtained with single-phase vapour-to-vapour, counter flow operation. First, the heat loss to the surroundings was quantified 
and found to affect the overall performance significantly. Second, the results from varying flow rates indicate that lower flow rates lead to higher 
outlet dew point values of the receiver stream which can be related to longer residence times. It was also found that moisture transfer is more 
strongly influenced by the flow rate through the receiver side than the donor side. Finally, five-point dew point profiles for both donor and receiver 
sides are reported for various temperature conditions. No stream wise variation in moisture flux was observed, and the average flux value increased 
from 3.3 x 107% kg s7! m7? at 30°C to 2.0 x 1074 kg s~! m~? at 70°C under fully humidified donor-side inlet conditions. 


© 2007 Elsevier B.V. All rights reserved. 
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1. Introduction 


One of the largest obstacles in the way of proton exchange 
membrane fuel cells (PEMFC) achieving commercial viability 
is the cost and size of the system. At least two ways to reduce cost 
and size are increasing power density (via water management, 
catalysts, materials, etc.) and trimming balance of plant costs. 
Reactant humidification subsystems are among the most expen- 
sive components in the balance of plant, and in addition can 
be a key performance enhancer. Technological improvements 
to reactant humidification will have a beneficial impact on the 
system power density and cost. 

Traditional stack humidification techniques such as satura- 
tion bubblers, direct injection, spray injection, etc., are being 


Abbreviations: DPAT, dew point approach temperature (°C) or (K); EE, 
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sible heat transfer effectiveness; WRR, water recovery ratio. 
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replaced by smaller, simpler, and more cost-effective solutions. 
Current stack humidification strategies are enthalpy wheels, self- 
humidification, and membrane humidifiers. Enthalpy wheels use 
a motor to rotate a desiccant-coated porous cylindrical core 
between the wet stream and dry stream. The core absorbs heat 
and moisture while exposed to the wet stream, and then cools 
and desorbs the moisture once rotated to the dry stream. This 
unit is advantageous because it is compact, effective, and has a 
very low pressure drop, but it is disadvantageous because it has 
numerous moving parts and seals, is expensive and demands an 
extra parasitic load to operate the rotary motor. 
Self-humidification methods [1,2] attempt to deliver water 
to the membrane using water or water-producing mechanisms 
available internally. Watanabe and co-workers have reported on 
self-humidified PEMFCs using a very thin membranes impreg- 
nated with particles of SiO2, TiO2, and Pt [1]. The oxides are 
highly hygroscopic and increase water retention. The noble 
catalyst enhances water production from Hag) and Og) that 
diffuse through the thin membrane and react internally. These 
authors maintain that the parasitic fuel losses are justified by 
the improved performance of the ionic conductor. However, 
the effects on membrane longevity and the membrane electrode 
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Nomenclature 

A membrane surface area (m?) 

Cp specific heat at constant pressure (J mol~! K7!) 

h enthalpy (Jkg~!) 

Jmem,H,0 average water flux across membrane 
(kg s7! m7?) 

L length of channel, 0.2 m (m) 

m mass flow rate (kg s7!) 

MW molecular weight (kg mol~!) 

Piso saturation vapour pressure (Pa) 

p total pressure (101,325 Pa unless specified other- 
wise) (Pa) 

pmo water vapour partial pressure (Pa) 

Q volumetric flow rate (SLPM) 

q' heat transfer rate to surroundings per unit length 
(Wm!) 

R? coefficient of determination 

T temperature (K or °C) 

Ta dew point above liquid water (K) 

tmem membrane thickness (m) 

U overall heat transfer coefficient (W m~? K7!) 

w effective perimeter for heat transfer to surround- 
ings (m) 

x distance along channel measured from dry-side 
inlet (wet-side inlet is at x= L) (m) 

Greek letters 

E effectiveness [0,1] 

Ymem membrane conductivity (W m`! K7!) 

À stoichiometric ratio 

E relative humidity [0,1] 

w humidity ratio, or specific humidity 
(kgH,0 Keay air) 

Subscripts 

H20 water vapour 

air air 

WI wet-side channel inlet (also used as an abbrevia- 
tion) 

WO wet-side channel outlet (also used as an abbrevi- 
ation) 

DI dry-side channel inlet (also used as an abbrevia- 
tion) 

DO dry-side channel outlet (also used as an abbrevi- 
ation) 

surr surroundings 

wet wet-side channel 

dry Dry-side channel 

x all probe locations along channel 


assembly response to sudden changes in current density are not 
clear. 

Although several conditioning schemes have been proposed 
to run PEMFC stacks on dry reactants, these schemes are only 


effective for low power applications (<5 kW) working at low 
temperatures (<60°C) [3]. Larger power applications require 
direct reactant humidification almost without exception. The 
reasons behind this requirement are not restricted to opera- 
tional constraints, and they can be explained by the fundamental 
properties of humidified gas mixtures. Larminie and Dicks con- 
sidered a PEMFC operating on dry reactants [3]. They treated 
the water vapour and the exiting streams as perfect gases, and 
assumed that all the product water was evaporated. With these 
assumptions, they calculated the partial pressure of water in 
terms of the total pressure at the cathode outlets, and the flow sto- 
ichiometry of dry air. Mérida reviewed these calculations [4] and 
used the resulting expressions to generate the curves in Fig. 1. 
These calculations illustrate that the stringent water require- 
ments within a PEMFC restrict operation to a very narrow range 
(i.e., a range for which & = 1 + ô, where d€ is small). Operating 
temperatures higher than 60 °C (which are desirable to minimise 
activation losses) correspond to very drying conditions for all 
practical flows at low pressures. 

A third type of humidification subsystem is a membrane 
humidifier. These passive devices recover humidity from the 
cathode exhaust and transfer it through a hygroscopic mem- 
brane to the cathode inlet stream (dry air). Two architectures 
are currently in use, both derived from compact heat exchanger 
designs. In the shell and tube design, the dry stream flows inside 
a bank of small membrane tubes while the wet stream flows 
over and around the bank of tubes. The plate and frame archi- 
tecture consists of membranes stacked on top of each other with 
a flow field plate separating the layers to allow flow. Membrane 
humidifiers are capable of a high moisture transfer capability 
at a reasonable pressure drop. This work characterizes a single 
channel version of a membrane humidifier. 

While humidifiers can be constructed with other membranes, 
Nafion membranes are the standard heat and mass transfer 
medium in a commercial membrane humidifier. The under- 
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Fig. 1. The variation of relative humidity at the cathode outlets as a function of 
temperature, stoichiometry, and operating pressure [4], adapted from [3]. 
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standing of water transport through Nafion has been the focus 
of intense analytical and experimental efforts. The work by 
Springer et al. in 1991 [5] was the first of many papers to model 
the diffusion of water as a phenomenon driven by a chemi- 
cal activity (water content) gradient. The water content of the 
membrane is calculated as a function of relative humidity from 
equilibrium sorption curves [6,7]. In practical fuel cell opera- 
tion, however, the membrane usually does not achieve sorption 
equilibrium, which can take between 100s and 1000s [8]. Only 
a few models consider the significance of surface mass transfer 
effects [8—11]. These models use a mass transfer conductance 
to relate the bulk gas water concentration to the water content 
on the membrane surface. The work of Ge et al. [8] goes a step 
further and suggests the absorption and desorption coefficients 
have a linear dependence on water content. 

The specific application of cathode humidifiers has not been 
given much attention in the literature until recently. In a notable 
early work Nguyen and White [12] discussed the effects of 
different humidification strategies on PEMFC performance. 
However, most of these techniques are outdated now. More 
recently, Chen and Peng [10] reported on a transient model of 
a planar humidifier in Simulink in the context of control sys- 
tems. Other relevant literature can be found in the field of energy 
recovery devices used to recuperate heat and humidity in build- 
ing ventilation systems [13-15]. Zhang and collaborators have 
extended heat transfer NTU analysis to mass transfer across 
membranes [14]. 

Experimental efforts to examine PEMFC membrane humid- 
ifiers have been sparse. Most experimental efforts involve 
condensing the outlet stream to obtain a time and space-averaged 
flux value for moisture transport across the membrane [16]. 
These flux values have been tabulated against geometric param- 
eters such as the ratio of residence time (e.g., the residence time 
that water molecules spend in the flow channel) to diffusion time 
(e.g., the time required to diffuse through air over the depth of 
the channel) [16,17]. This ratio has been used successfully to 
specify the optimal operating ranges for planar humidifiers and 
to convert the windows in parameter space into design specifi- 
cations. Other experimental work was provided by Park et al. 
in 2005 [18], who reported data for liquid water-to-gas internal 
humidification using multiple membranes. 

In this work, we report the effects of flow rates on perfor- 
mance of planar membrane humidifiers operating under single 
phase conditions. We then address a gap between simulations 
based on analytical models and the macroscopic experimen- 
tal studies. The former provide complete spatial mapping of 
humidity and temperature while the latter rely on time and 
space-averaged fluxes. We provide humidity measurements at 
different locations along the channels to assist ongoing modeling 
efforts. 


2. Experimental 
2.1. Setup 


The experimental setup consisted of a single channel mem- 
brane humidifier as illustrated in Fig. 2. A 1 cm x 1 cm x 20cm 
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Fig. 2. Single channel humidifier. 


flow channel was machined into two plates, which were then 
mated together with Nafion 117 membrane as the hygroscopic 
separator. The Nafion was pretreated by boiling in a 3 wt.% 
H20, solution for 2h, rinsing with de-ionized water, boiling 
in 1.0M H2S0; for 2h, and finally rinsing and storing in de- 
ionized water. The membrane was removed from de-ionized 
water and exposed to ambient room temperature and humid- 
ity conditions for at least 24h prior to testing. Five evenly 
spaced 1/16in. probe holes were drilled into the flow chan- 
nels for thermocouple or dew point sensor placement. Calibrated 
type-T unsheathed thermocouples were used due to their reason- 
able accuracy, suitability in moist air environments, and small 
size (to minimize interference with the flow distribution in the 
channel). 

A Vaisala HMT337 dew point transmitter was used to 
measure moisture content. Because this is a polymer-based 
capacitive sensor, an effort was made to ensure the probe tem- 
perature was equal to the process temperature. The dew point 
probe was placed in a cavity, and at the time of measurement 
a sample of gas from the channel was allowed to flow through 
this cavity to perform accurate dew point measurements. The 
Vaisala meter measures relative humidity at +2% and temper- 
ature at +0.2°C. The Hyland and Wexler equation was used 
throughout this work to relate saturated water vapour pressure to 
temperature: 


In pifo = a + C2 +CT + C4T? + CsT? + C6ln(T) (1) 
where T is in Kelvin and Pio is in Pa. The constant coefficients 
are listed in Table 1. 

In this work, dew points below 0°C refer to the temper- 
ature at which the parcel of moist air would be saturated 
above super cooled liquid water. This is distinct from the 
frost point, which correlates moisture content to the satura- 
tion pressure above a frozen surface (i.e., ice). In all cases, 


Table 1 

Coefficients in Hyland and Wexler saturation vapour pressure equation 
C1 —0.58002206 x 104 

C2 0.13914993 x 10! 

c3 —0.48640239 x 107! 
C4 0.41764768 x 1074 

C5 —0.14452093 x 1077 
C6 0.65459673 x 10! 
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Fig. 3. Experimental setup schematic. 


however, negative dew point values indicate very low water 
concentrations. 

Fig. 3 shows the experimental setup. The dry side stream (or 
receiver side) was delivered through a mass flow controller and 
heated to the desired inlet temperature, Tpr. The wet side stream 
(or donor side) passed through a mass flow controller, a satu- 
ration point bubbler humidifier, and then heated to the desired 
temperature, Twr. The moisture content of this stream was con- 
trolled by the temperature and pressure in the bubbler. In actual 
operation, the wet side stream, coming from the cathode exhaust, 
would be oxygen-depleted moist air. In this work, regular air was 
used on both sides. 


2.2. Methodology 


Three related questions were addressed in this work. First, in 
our configuration, the heat lost to the surroundings was compa- 
rable to the enthalpic flow through the humidifier and thus not 
negligible. The first step then was to characterize and quantify 
the heat loss of the humidifier. A stream may gain or lose heat 
energy through the channel walls to the surroundings, through 
the membrane to the opposite stream, and through latent heat 
transfers in the presence of phase changes. Single phase condi- 
tions were imposed throughout this work for simplicity, although 
in practical operating conditions there may be condensation on 
the wet side. 

Next, the effects of flow rates were studied. The effect of 
flow rates on humidifier performance is of utmost interest in 
humidifier design because it characterizes the humidifier just 
as a polarization curve characterizes a fuel cell. The fuel cell 
system’s air flow rate will be determined by the stoichiomet- 
ric ratio and current demands, so it is desirable to characterize 
the humidifier’s moisture-transferring ability as a function of 
flow rate. Furthermore, membrane humidifiers are currently 
passive devices, humidifying the reactant as much as pos- 
sible before entering the cell. However, in some operating 
regimes or future applications, it may be desirable to control 
the humidity over a range. One example of a control mecha- 
nism is a dry air bypass gate as suggested in [10]. In this case, 
there would be a differential flow on either side of the mem- 


brane. It would then be advantageous to fully understand the 
tradeoffs of changing wet side and dry side flow rates individu- 
ally. 

A randomized full-factorial experiment was carried out with 
three replications at three levels of both wet side and dry side flow 
rates under four sets of conditions. The experimental conditions 
for each case are summarized in Table 2. 

The third portion of this study aimed to provide information 
regarding the stream-wise variation in moisture flux. Knowing 
where most of the moisture transfer occurs in counter flow mode 
could provide valuable information leading to optimization of 
the channel length. Two replications with randomization were 
taken at each of the 10 probe locations. This was performed for 
the same four sets of conditions. 

The seven independent variables that can be varied in this 
arrangement are: dry inlet temperature, dew point, and flow 
rate (Tpr, Tap1, and Qairpr), wet side inlet temperature, dew 
point, and flow rate (Twr, Tawi. and Qair wi) and the temper- 
ature of the surroundings (Tsurr). The dependent variables can 
include temperatures or dew points at the 10 probe locations, 
Ty or Tg... Due to the high number of variables, four individual 
cases representing different operating regimes were considered. 
The case definitions and experimental testing approach is shown 
in Table 2: Testing Matrix. 


2.3. Measurement technique 


Fig. 4 illustrates the dew point measurement technique. In 
steady state, if the convective velocity of the air/water mixture 
along the channel is greater than the diffusive velocity of water 
in air across the channel, then a significant concentration gradi- 
ent in the y-direction will develop. As a result, the sample probe 
location on the wall opposite the membrane may measure an 
artificially low concentration value on the dry side and an artifi- 
cially high concentration value on the wet side. This effect will 
become more noticeable further along the channel were con- 
centration gradients grow larger and the mixture is no longer 
well-mixed. To mitigate this effect, the channel exit stream was 
partially restricted during measurements to force a larger sample 
of air/water mixture through the sensing cavity. 
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Table 2 
Testing matrix 


Case Research focus 


Fixed temperatures in 
(°C) flows in (SLPM) 


Independent Dependent 


Case 1: non-isothermal case Heat loss to surroundings 


Effect of flow rates on moisture transfer 


Streamwise variation in moisture flux 


Effect of flow rates on moisture transfer 
under isothermal conditions 


Cases 2, 3, and 4°: isothermal 
cases at low, medium, and 
high temperatures, 
respectively 


Streamwise variation in moisture flux 


Tp = 85 
Tap = —40 
Qair, pr = 1.0 
Tw = 85 
Tawi = 25° 
Qair,wi = 1.0 
Turr =21 
Tp = 25 
Tap =-40 
Twr = 80 
Tawi = 25° 
Tsurr = 21 
As above, with 
Qair.p1 = 1.0 
Qair,wi = 1.0 


Tot = Tsurr 


Streamwise location Ty 


0.4< Qair,DI <1.0 
0.4< Qair,WI <1.0 


Ta,DO 


Streamwise location Tx,Tax 


0.4< Qair, DI < 1.0 Tapo 


Tapi = —40 
Twi = Turr 
Tawi = Tsurr 


0.4< Qair, WI < 1.0 


Where 

Tsurr = 30 (Case 2) 
Tsurr = 50 (Case 3) 
Tsurr = 70 (Case 4) 
As above, with 
Qair,pi = 1.0 
Qair,w1 = 1.0 


Streamwise location Ta x 


à This is the maximum dew point allowable to ensure no condensation occurs in the wet side as the temperature drops; limited by humidifier heat loss to the 


surroundings. 
P These cases use an oven to impose isothermal conditions. 


3. Theory 
3.1. Flow correction 


As a result of the experimental set up shown in Fig. 3, one 
must consider the addition of water when stating the actual wet 
side flow rate delivered to the membrane humidifier. The amount 
of added water is quantified, assuming an ideal gas mixture and 


Sensor reads 
artificially low on 
dry side and high 

on wet side 


(Dry Side Inlet / 
Wet Side Outlet) 


Dry Side Outlet) 


Fig. 4. The measurement technique showing conceptually the concentration 
gradients across the channel and their effect on measurement precision. 


100% effective bubbler, via: 


sat 
z Pho MWmn,o 
H0 = 
p — pio MWair 


Mair (2) 


where Pito is the saturation pressure at the bubbler tempera- 


ture, and p is the total pressure in the bubbler. The mass flow 
rate of water and air are summed to determine the actual flow 
composition and rate delivered to the membrane humidifier. For 
reference, Fig. 5 has been included to show the appropriate cor- 
rection factors (either volumetric or gravimetric) to obtain the 
total flow rate of the mixture. Higher backpressure build up at 
the saturation bubbler results in less water being evaporated by 
the flow. Should one not account for these corrections, the max- 
imum error will be found in cases with high dew point and high 
flow rates (which cause higher back pressures). The sensitivity 
to dew point is much higher than to bubbler pressure. 


3.2. Performance metrics 


Several candidate metrics to measure a humidifier’s moisture- 
transferring ability have been used in the past. The stream of 
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Fig. 5. Volumetric and mass flow rate correction factors due to addition of water 
for various bubbler backpressures. 


interest is the dry side outlet which will be delivered to the cath- 
ode and the variable of interest is the amount of water carried in 
that stream. An appropriate performance metric often depends 
on whether one is comparing different humidifiers at identical 
conditions or comparing the same humidifier at different con- 
ditions. Generally, however, outlet relative humidity and outlet 
dew point are poor choices because the temperature (if using 
relative humidity), pressure, and flow rate are still needed to 
quantify the flow rate of water present in the stream. Further- 
more, there is no indication of the maximum value that could 
have been attained under the conditions. 

The dew point approach temperature (DPAT) is the differ- 
ence in dew point between wet side inlet and dry side outlet—a 
lower DPAT is indicative of better performance. When compar- 
ing different humidifiers or the same humidifier at conditions 
with constant wet side inlet dew point, the DPAT provides a 
practical evaluation of the performance. To exemplify an area 
of caution when using DPAT, consider the case where a non- 
perfect humidifier is run with Tg wr=40°C and yields a DPAT 
of 5 °C. If it is then run under the same flows with Ty wy = 80 °C 
and yields a DPAT of 30 °C, the DPAT performance is worse but 
it has in fact transferred more water. Also, DPAT is a skewed 
scale due to the nature of the saturation pressure curve; a 3 °C 
DPAT is not twice as good as a 6 °C DPAT. 

The total water transfer is calculated as the difference in 
amount of water between inlet and outlet on either stream. The 
average water flux is the total water transfer normalized by total 
membrane area: 


MH,0,DO — MH)O,DI 


A (3) 


Jmem, H0 = 


where the mass flow rates of water, 774,0,po and my,0,p1 are 
calculated using Eq. (2) with Piso replaced with the water 
vapour partial pressure, py,o. Either the total water transfer 
or average water flux are suitable metrics when comparing 
the same humidifier at different operating conditions, but only 
the total water transfer is suitable when comparing differ- 


ent humidifiers. Note, however, that these metrics have no 


indication of the maximum amount that could have been trans- 
ferred. 

There are several first Law effectiveness measures which nor- 
malize humidifier performance on a scale from 0 to 1. The WRR 
is defined as the ratio of the total water transferred to the quantity 
of water available in the wet stream: 


MH,0,DO — MH)O,DI 
WRR = 2 2 


(4) 
™¥,0,WI 

The WRR is only suitable if running on dry reactants, other- 

wise the maximum amount that could be transferred is less than 

mm o,wi. The sensible, latent (or moisture), and enthalpy (or 

total) effectiveness factors for heat and moisture transfer devices 

are defined as: 


_ Mair,p1(Xpo — Xnr) 
(hair)min(Xwi — Xp1) 


(5) 


where X is either the temperature (7), humidity ratio (w), or 
enthalpy (h) depending on whether the sensible effectiveness 
(SE), latent effectiveness (LE), or enthalpy effectiveness (EE) is 
being calculated. The SE is a measure of heat transfer; not mass 
transfer, and therefore is of limited use in humidifier design. 
The LE and EE are useful measures and correct the WRR in the 
case of partially humidified dry side inlet conditions. EE can be 
misleading from a moisture transfer point of view if the sensible 
heat exchange makes up the bulk of the total enthalpy exchange. 
The LE is the better of the three to isolate moisture-transferring 
effectiveness. 

Finally, a deficiency of all of these effectiveness parameters 
is that they only consider mass and energy conservation. Second 
law considerations such as irreversibilities due to heat transfer, 
phase changes, absorption/desorption, and friction losses are not 
accounted for so a thermodynamic maximum performance is 
difficult to predict. 


4. Results and discussion 
4.1. Heat loss to surroundings 


The amount of heat lost to the surroundings was quantified 
by how it affects the temperature profiles in co-flow operation 
with both inlet streams flowing dry air at the same temperature. 
Neglecting changes in potential or kinetic energy, the heat loss 
per unit length for this configuration, q’ (W m7 !), was calculated 
according to the equation: 

i d(mcpT) _ dT 


= = h 6 
q dx Map dx (6) 


The derivative of the quadratic fit (R? = 1.00) shown in Fig. 6, 
the appropriate flow rate, and the specific heat capacity of dry air 
(cp =1007J kg"! KT!) were used to calculate the heat loss rate. 
If the heat loss to the surroundings is modeled with an overall 
heat transfer coefficient, the following equation applies [19]: 


q = Uw(Tsur — T) (7) 


And the product of overall heat transfer coefficient, U, and 
effective perimeter, w, can be extracted from the slope of a 
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Fig. 6. The temperature profile and the heat loss rate to surroundings with dry 
air in coflow mode (1.0 SLPM). 


linear fit in a plot of g’ versus (Tsun — T) (Fig. 7). A value of 
Uw = 0.149 + 0.023 WK! m`! at a 95% confidence inter- 
val was obtained for this particular experimental setup. The 
Uw term can be assumed constant for fully developed laminar 
flows because the Nusselt numbers, thermal conductivity of air, 
and consequently convective heat transfer coefficients all remain 
approximately constant. 


4.2. Effect of flow rates on humidification performance 


Reporting solely the outlet dew point or relative humidity 
[18] as the performance variable can be misleading because it 
is possible and common to find that a lower outlet dew point 
at a higher flow rate can represent more mass of water trans- 
ferred than a higher dew point reading at a lower flow rate. 
Therefore, the average water flux across the membrane and 
WRR performance indicators are also considered in the cur- 
rent section. The WRR is suitable as an effectiveness measure 
because the dry side inlet was dry air in all cases. Contours of 
these three metrics for the ranges of flow studied are shown in 
Figs. 8-11. 
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Fig. 7. Determination of overall heat transfer coefficient. 


Two observations can be made with regard to Figs. 8a, 9a, 
10a, and 1| 1a. First, in all cases the outlet dew point increases at 
lower flow rates. This result might be expected since lower flow 
rates give both the wet and dry gases a longer residence time 
in the channel, allowing for more moisture to be evaporated. 
Second, the effect of the dry side flow rate on outlet dew point 
is generally much more pronounced than the wet side flow rate 
(evidenced by the near-vertical contour lines). This suggests that 
an abundance of moisture exists at the membrane interface on 
the dry side and it is only the dry side’s time in the channel that 
limits how much water is evaporated. 

Figs. 8b, 9b, 10b, and 11b uncover no clear trends from case 
to case. However, it should be noted that a higher outlet dew 
point does not necessarily mean that more moisture has been 
transferred (in many cases the opposite is true). This is an impor- 
tant consideration when comparing performance data because 
using dew point alone may lead to erroneous interpretations. For 
example, Park et al. [18] concluded that flux across the mem- 
brane increases with an increase in flow rates, whereas Ge et 
al. [8], in a similar application, concluded that the outlet dew 
point decreases with increasing flow rates. The two conclusions 


Case 1 - Di: T=20°C, Ty= -40°C / WI: T=80°C, Ty=25°C / T,, =21°C 
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Fig. 8. Case 1 (non-isothermal)—effect of flow rates on (a) outlet dew point, (b) average water flux, and (c) water recovery ratio. 
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Case 2 - DI: T=30 °C, T= -40°C / WI: T=30°C, T,=30°C / Teun =30°C 
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Fig. 10. Case 3 (isothermal, 50 °C)—effect of flow rates on (a) outlet dew point, (b) average water flux, and (c) water recovery ratio. 


appear in conflict, but are solely the result of the metric used 
to indicate performance. A second observation with respect to 
these graphs is that there exists an optimal combination of flow 
rates that yields the highest water flux. However, this optimum 
varies with operating conditions in an unclear way and in some 
cases is not reliably discernible with the data presented here 
(Fig. 11b). Nonetheless, the importance of such information can 


Case 4 - DI: T=70°C, T= -40°C / WI: T=70°C, T,=70°C / 
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humidifier. 
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Fig. 11. Case 4 (isothermal, 70 °C)—effect of flow rates on (a) outlet dew point, (b) average water flux, and (c) water recovery ratio. 


416 


Table 3 
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Calculated accuracy and precision estimates for dew point readings in factorial experiments 


Wet side (SLPM) 


Accuracy? (°C) dry side flow rate (SLPM) 


Precision? (+°C) dry side flow rate (SLPM) 


0.4 0.7 1.0 0.4 0.7 1.0 

Case 1 

0.4 1.38 1.66 1.94 1.16 1.37 0.50 

0.7 1.41 1.65 2.07 0.75 0.79 0.05 

1.0 1.47 1.79 2.12 0.02 0.86 0.18 
Case 2 

0.4 0.80 1.10 1.20 2.07 2.12 2.89 

0.7 0.71 0.99 1.18 1.89 1.73 1.45 

1.0 0.70 0.91 1.17 0.70 0.98 0.24 
Case 3 

0.4 0.90 1.32 1.44 0.61 0.95 3.16 

0.7 0.86 1.09 1.37 1.86 2.50 2.05 

1.0 0.85 1.00 1.33 0.69 1.75 1.14 
Case 4 

0.4 1.22 1.46 1.71 2.59 2.02 1.65 

0.7 1.10 1.51 1.70 0.66 3.08 0.98 

1.0 1.20 1.49 1:72 0.84 2.54 1.75 


* Based on manufacturer’s published accuracy limits for given conditions. 


P Based on standard deviation of three replications. 


All cases exhibit similar trends on the water recovery ratio 
graphs; better effectiveness at high dry side flow rates and 
low wet side flow rates. Note, however, that a high effective- 
ness rating (a relative measure) does not necessarily mean the 
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Fig. 12. Moisture profiles—Case | (non-isothermal). 
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Fig. 14. Moisture profiles—Case 3 (isothermal, 50°C). 
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70°. The reason for this is because the higher temperature dew With the net change in water flow rate presented in the mar- 


points carry a lot more water into the system, but the humidifier gins of Figs.12b—15b, it is possible to compare the amount of 
is sized too small to transfer a significant portion of it to the dry water lost from the wet side to the amount of water gained 
side. by the dry side. Although the difference was only statistically 

The estimated accuracy (based on sensor accuracy) and pre- significant in Case 1, the amount of water ‘lost’ by the wet 
cision (standard deviation between replications) of these results side is slightly higher than the water being ‘gained’ by the 
are summarized in Table 3. At higher temperatures, even small dry side stream. The possibility of a leak is unlikely; a more 
error in measurement can result in large errors in the calculated probable cause may be a result of measurement technique. As 
values of average water flux and water recovery ratio. 


Table 4 
4.3. Stream wise variation in moisture flux Calculated average water fluxes from Figs. 12b, 13b, 14b, 15b 
Calculated flux value Statistically 

The measured dew point profiles and calculated mass flow of (x10~> kg s~! m~?) different slopes? 
water (Eq. (2) with Piso replaced with the water vapour partial Case 1 
pressure py,0) are shown in Figs. 12—15. For these figures, the Wet side 0.51 Yes 
arrows indicate direction of flow and the error bars shown are the Dry side 0.26 
manufacturer’s stated accuracy limits or the standard deviation Case 2 
of the replications, whichever is higher. Wet side 3.60 No 

In all cases, the mass flow rate of water in the channels Dry side 3.04 
(Figs. 12b, 13b, 14b, and 15b) appears to change linearly Case 3 
along the channel. Table 4 summarizes the fitted linear models. Wet side 8.25 No 
Although no phenomenological significance can be ascribed to Dry side 7.76 
the linear profile, it provides a useful empirical correlation for Case 4 
engineering design. A linear profile implies that the flux across Wet side ce No 

Dry side 18. 


the membrane is approximately constant along the length of the 
channel. è At a 0.05 significance level, determined by t-test. 
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mentioned in Section 2.3, the effect of concentration gradient 
across the channel may cause dry side measurements to be 
lower and wet side measurements to be higher than their true 
values. 


5. Conclusions 


An experimental investigation into the effect of operating 
conditions on a single channel gas-to-gas Nafion membrane 
humidifier has been carried out. A method for quantify- 
ing the heat loss rate to the surroundings per unit length 
was demonstrated by calculating an overall sensible heat 
transfer coefficient. For this particular geometry, a value of 
0.149 W K7! m7! was determined. This value is expected to 
change slightly with humid air, and change markedly in the 
presence of condensation on the walls. 

The effect of flow rates on humidifier performance was stud- 
ied by varying the dry side and wet side flow rates individually. 
In reporting results pertaining to moisture transfer, trends will 
vary depending on the metric used to indicate performance. The 
outlet dew point increases with lower dry side flow rates but 
was not significantly affected by wet side flow rates. The aver- 
age water flux, however, did not follow the same trends since 
a higher outlet dew point or relative humidity seldom trans- 
lates into better moisture transfer when comparing different flow 
rates. Certain flow rate combinations yielded optimum water 
flux across the membrane at each temperature, but these opti- 
mum flow rate combinations appear to be condition-dependent 
because no clear trend was uncovered with increasing tempera- 
ture. 

Finally, stream-wise dew point and mass flow of water pro- 
files were presented. In the cases studied, moisture flux across the 
membrane did not vary significantly in the stream wise direction, 
and increased with higher operating temperatures. 
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